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Abstract 



The article deals with the subject of solving the problem of canonical- map syn- 
qq ■ thesis for Hamiltonian systems. For this purpose, the controlling-function method 

^sO ■ has been developed that allows appropriate changes of the variables in terms of 

£^ | calculus of canonical variation, starting from their target conditions. To use the 

canonical formalism, the initial dynamic system that employs changing Lagrange 
■ multipliers is reduced to a Hamiltonian system in an expanded phase space, fol- 

. lowed by the construction of controlling function. The algorithm suggested for the 

canonization of controlled mappings has an advantage over the known procedures, 
and first of all, redundancy in the procedure that chooses regulated coordinate 
^ , transformations as a base for a goal-seeking synthesis scheme. 



Introduction 



Below we suggest a new procedure of variational changes in dynamical and mechanical sys- 
tems, which is named the method of controlled (or synthesized) mappings. In distinction 
to the known methods of generative functions, those of Lie generators, parametrization 
in the Hilbert-Courant form, etc. [1-5], the method allows a large extent of general- 
ity and versatility (for example, with respect to the transformations of coordinate axes, 
etc.) because it makes possible to synthesize mappings in dependence on the chosen tar- 
get conditions without resort to some laborious-intensive procedures (as it takes place in 
a classical variant of the generative- function method). 

The formalism of generative functions unfortunately does not permit beforehand a new 
system of variables suitable to finding solution of the selected problems (normalization, 
integration, and so on). On the contrary, the method of controlled mapping makes possible 
to solve these problems constructively basing on the beforehand formulated conditions. 
It might be caused by some shortcomings in the method itself, which presented, as so 
often was the case, little more than derivation of the related equations in the required 
quantities, but did not present their analytical solution. 

In Section 1, a controlled change of variables in dynamic systems has been considered. 
For this purpose, canonical equations are derived variationally, with giving proof of a nec- 
essary Weierstrass condition for a minimal action functional in terms of invariant Hilbert 
integral. The controlled mapping is synthesized in phase space by means of a controlling 
function, and its discovery becomes the most important problem of the method. 

In Section 2 the most distinctive singularities of the suggested controlled transforma- 
tions are pointed. Here we, first of all, reveal the conditions for appearing the properties 
of canonicity and invariance. Various forms of canonicity criteria for the controlled map- 
pings are found, and some possible ways of their support are discussed. Great attention is 
paid to the derivation of Hamilton- Jacobi equation for the controlling potential. For the 
introduced controlled mappings, the invariance properties are established. 

Section 3 presents description of energy transformations in an initial dynamic system, 
which were induced by the canonical controlled mappings. The relationship between the 
acting and controlling functions is established, which provides the controlling transforma- 
tion field in phase space. By a basic target condition imposed on the new Hamiltonian, 
they present an algorithmic scheme for finding the controlling function. 

Section 4 is small in volume but important in theoretical generalizations. The direct 
analogy between small controlled canonical mappings and infinitesimal canonical trans- 
formations is firstly revealed. It is secondly shown that the canonization of controlled 
transformations provides target conditions without using a special choice of initial data. 

1 Controlled change of variables 

The controlled transformation of variables in any moving dynamic system is presumed 
to be a change of system variables in the process of system movement as time passes, i. e. 
supplying a controlled signal in the form of controlled mapping to an input converter of 
initial variables x(t). 



Let movement of the dynamic system be described by an n-dimensional vector equation 



x = f(x,t), xeR n , (1.1) 

where x = (x±, x n ) is the phase vector, Xi = Xi(t), i — l,n are phase variables (coordi- 
nates, velocities), t G [to, ti] C R is time. It is thought that /, df/dx, df/dt e C 1 [to ) ti ] 
on the right side of equation (1.1) satisfy the existence and uniqueness condition applied 
to the Cauchy problem in some limited region D C R n x i? , i. e. equation (1.1) is satisfied 
by a unique continuous integral path in the expanded phase space R n x R that originates 
at the point x^ = Xi(to) and terminates at the point Xu = Xj(ti), % — l,n. 

Given differential constraints in the form of equation of motion (1.1), let us introduce 
the Lagrange multipliers Aj(t), % — l,n, where \(t) is the unknown temporal functions 
determined as part of the solution of optimization problem for the action functional S, 
with the provision of constrained equation (1.1). The Lagrange multipliers can be seen 
as weighting action-functional integral multipliers 

rti 

S= / L(x,x,X,t)dt, (1.2) 

Jto 

where L — X (x — f) is the Lagrange function (Lagrangian). 

The variation of variable quantity is, as we know, the difference of values of this 
quantity on a permissible reference trajectory and on an appropriate permissible trajectory 
called the comparison trajectory. The variation of the quantity calculated for the same 
moment of time t is called the isochronous variation and denoted by 5. If in calculating 
the variation the value of the quantity is considered on the comparison trajectory at time 
t + At, this variation is called a total variation. Let the total variation be denoted by A. 

Supply also the quantities taken on the comparison trajectory by a line over. Then 
the relationship between the total and isochronous variations is specified by the equalities 

Axi = Xi(t + At) - Xi{t) = Xi(t + At) - Xi(t) + Sx u 



Axi = 5xi + Xi At, i = l,n, (1.3) 

and the quantity Xi in formula (1.3) substitutes for Xi if they were considered to be 
diverged infinitesimally. It is here thought that we have the permutation relations Sii = 
Xi — Xi = d (xi — Xi) / dt — d (5xi) / dt. 

To seek AS, let us go from the Lagrange function L to the corresponding Hamilton 
function H using the Legendre transformation 

dL 

H(x, X,t) — —x — L, L — L(x, x, A, t). (1.4) 
ox 

According to the definition of function L and by (1.4), we obtain 

dL 

8i =X - » = Xi - (L5) 



The vector multiplication is here and further on considered as a scalar product. 



With relations (1.4), (1.5) we can write [6-8] 

rti+Aii _ 



dt 



AS = Xxdt- Xxdt- A Hi 

Jto+Ato J t J t 

= / (Xx — Xx) dt + Xx \ ti Ati — Xx \ t At 

J to 

- f h SHdt-^Ah + H At , (1.6) 

Jto 

where H = H \ t=to , H 1 — H \ t=tl , A — A = 5X, and all equalities in (1.6) are valid accurate 
within infinitesimals of the order of smallness higher than first. The first summand in the 
right side of formula (1.6) may be written with the mentioned accuracy in the form 

rti _ _ rti _ rti 

/ (Xx — Xx) dt — X(x — x) dt + xdXdt 

Jto Jto Jto 

= [ h X — (5x)dt+ f h x5Xdt= f t \-X5x + x5X)dt + XSxf' (1.7) 

Jto dt Jto Jto 

Using relations (1.3), (1.7), and 

5H(x, X,t) = ^5x + ^-5X 
ox oX 

for the calculation of AS being accurate with the second and higher orders of smallness, 
we obtain 



AS = - 



^\ fdH -\ . fdH 
^ fe+ 



to 



dt 



+ (XAx-HAt)\l\ (1.8) 

If there is an extremal of the functional S, it is necessary to fulfil the stationary 
condition AS* = 0. When the variations in the integral (1.8) are independent and when 
we take account to the equalities (1.5), we obtain necessary conditions in the form of: 

(1) Euler equations (Euler-Lagrange) with respect to the multipliers A 

• dH " df k . 

>H = -dx- = ~h Xk dx~' * = 1 ' n; (L9) 

(2) equations of motion (1.1) with respect to the variables x in the form of Euler equa- 
tions 

dH 



Xi = — = fi i = l,n; (1.10) 



(3) the universality condition 



(XAx-HAt)\l=0. (1.11) 



I to 



It is more believed that the right sides of equations (1.9), (1.10) satisfy the existence 
and unique conditions for the Cauchy problem in a limited region D, namely, they satisfy 
the continuous vector-function f(x,t) and limited elements of the Jacobi matrix f x (x,t). 
These equations therefore fit a unique continuous integral trajectory in the expanded 
space R n x R that issues out of the point with coordinates xq, Ao at the initial moment 
of time t . 

Note that Euler-Lagrange equations (1.9), A = — f x X, f x = df/dx, are the Euler 
equations of variational Lagrange problem (1.2) with the fixed ends that serve to finding 
an extremal of the functional S, 

L — —L- 

1, XI 

at 

where ^ ^ 

L x — ~t\ = — fx\ L x = -— = A. 

ox ox 

It is clear that the system of equations (1.9), (1.10) is well-evident to form a canonical 
system with the Hamiltonian H, coordinates Xi, and factors Aj. This system is different 
from a familiar canonical system by no more than the linear dependence of H on A j . What 
is more, for a total derivative of the function H = H(x, A, t) with respect to the time, we, 
according to (1.9) and (1.10), have 

dH dH dH ^ dH 

dt dx dX dt 

dH dH dH dH dH _ dH 

dx dX dX dx dt dt 



Because by relations (1.5) the expression 

^ = Xf + Xf = Xf + Xf x x + Xf t = Xf t 

is valid, we shall from here find that on the extremal obeying the necessary conditions for 
the extremal the equality f x = df/dx, f t = df /dt takes place, we from here find that, 
the equality 

r ti 

H = h+ Xft dt,h = const (1.12) 



to 

must be fulfilled. If the vector function / does not explicitly depend on t, it follows 
evidently from (1.12) that H = h. 

We pointed out that from equations (1.9), (1.10) their contingency follows, where the 
vector X(t) was that of conjugate variables. But other than these, the functions L and H 
are conjugate characteristic functions derived from the Legendre transformations because, 
by (1.4), (1.5), and (1.10) on the extremal trajectories of functional (1.2) while differential 
relation (1.1) is provided, the following equalities are resulting: 

flr f)M 

L = —x-H = X — -H = Xf-Xf = 0. 
dx dX 



We also further on need certain results of variational analysis concerning the notion of 
integral invariance. The calculus of variation [9] for problem (1.2) with the Lagrangian L 



the class of curves C with the ends given on the interval [to, t±\ is considered. It is necessary 
to find a minimum of the integral S(C) (1.2) along C, where x(t) is the representation 
of curve C. The integral is here taken over the interval of time [£o,£i], within which the 
curve x{t) is defined. 

It is stated that if M is the class of curves with the given ends A = x(t ), B = x(t\), 
and C is a curve from this class, the quantity S(C ) is their minimum S(C) with respect 
to C G M on the assumption that there exists a total derivative $ of the function 
= ty(x, A, t) such that L = $ along the curve Co and L > $ along all other curves from 
the class M. 

Explain that we here have $ = \I/ = \P X x + ^\ A + \I/ t . The statement in question may 
be proved straightforwardly: VC G M the value S(C) is as follows: 



S(C) > [ <5>dt = I Ldt = S(C ) 

J(C) J (Co) 



1(C) J (Co) 

Definition 1 The vector function g = g(x, t) is called a geodesic inclination if there exists 
a total derivative $ = $ = + $aA + $t of the function ^ = ty(x,\,t) such that 
V (x, x, A, t), where (x,t) are lying in the region of definition of the function g(x,t), the 
inequality 

L > $ 

is valid, and L = $ at x — g. Here we have g(x, t) : R n x R — » R n . 

Thus, min (L — $) = is achieved at x — g. Providing the coupling equation (1.1), we 
therefore have 

L = A (x - f) = S> x x + * A A + * t> 
and we from where obtain a chain of equalities 

V X = L ± , q x = L x = 0, % = L-L ± x = -H. (1.13) 

Let us add to definition 1 [9] that if g is the geodesic inclination, the geodesic inclination 
curves are called solutions of the differential equations x — g, where g = g(x, t). The family 
of such curves is called a geodesic family. The geodesic family is named the geodesic flow 
if it covers the region, where the function g is only once defined. 

Note that in the agreed notation, the integral 

S=f [L(x,g,X,t)-(x-g)L g (x,g,X,t)]dt (1.14) 

J(C) 

named the invariant Hilbert integral, where L g = dL/dg, is independent of the integration 
path C, but only dependent on the end of the curve C. Due to relations (1.13) at x — g, 
integral (1.14) is indeed equal to the value f,™ $ dt = ^f(B) — ^(A), where A and B are 
the beginning end of the integration path; in our case $ dt = ^f(B) — ^f(A). 

With the use of inequality of the form L — $ > and the substitution of $ as 
a integrand (1.14), we obtain the Weierstrass condition E > 0, where the Weierstrass 
function E = E(x, g, x, A, t) is of the form 



E = L(x, x, A, t) - L(x, g, A, *)-(£- g) L g (x, g, A, t). 



(1.15) 



The min (L — <3>) = is thus achieved at x — g. When relation equation (1.1) is 
satisfied, we therefore have 



L = A (x - f) = ^ x x + ^ A A + %, 

we from where obtain the chain of equalities 

V X = L ± , q x = L x = 0, ^ t = L-L ± x = -H. (1.13) 

Add to the definition 1 [9] that if x — g is the geodesic inclination, the solutions of 
the geodesic equation x = 0, where g = g(x,t), are referred to as the geodesic inclination 
curves. The family of such curves is called a geodesic family. The geodesic family is 
referred to as a geodesic flow if it only once covers the region of definition for the function 
9 = 9(x,t). 

Point out that in the agreed notation the integral 

S= [ [L{x,g,X,t)-{x-g)L g {x,g,X,t)\dt (1.14) 

J(C) 

named the Hilbert integral, where L g = dL/dg, is independent of the integration path 
C, but dependent on the ends of the curve C. By (1.13), the integral (1.13) at x — g is 
indeed equal to $ dt = ^(B) — ^(A) at x — g, where A and B are the initial and 
terminal ends of integral pathes, being L(x, g, A, t) — A (g — f) in our case. 

Using the inequality L — $ > from definition 1 and substituting the integrand 
(2.14) for <3>, we obtain the Weierstrass condition E > 0, where the Weierstrass function 
E = E(x, g, x, A, t) has the form 

E = L(x, x, A, t) - L(x, g, A, t) - (x - g) L g (x, g, A, t). (1.15) 

Point to the fact that if g is the geodesic inclination, we have L = $ + E, hence, if 
VC 6 M, it joints two levels ^ = and ^ = ^2, we from here obtain the Weierstrass 
formula of calculus of variation 



S(C) =* 2 -*i + / Edt, 

J(C) 



or 



S(C) = S(C )+ f Edt, 

J(C) 

where Co is the curve of flow, and C is any curve joining the levels of the curve Co- 

A necessary Weierstrass condition (E > 0) for the minimum of the action functional 
S (see (1.2)) can be probed as follows. When g = g(x,t) substitutes for the function x, 
the Lagrange function L(x,x, X,t) changes into the function L(x,g, X,t). Expanding into 
series, we find 



dL 

L(x, x, A, t) = L(x, g, X,t) + — 



(x-g) +a{\x-g\), 

x=g 



where a(-) denotes a set of infinitesimals of the order higher than first order with respect 
to | x — g |. The written expression gives the Weierstrass function E (1.15) if we neglect a. 



We have for the system (1.1) in the corresponding variational problem, according to 
formulas (1.4) and (1.5), 



dL 

L = \[x-f{x,t)\, ^- = A. 

Hence, 

E = X(g-f)-X(x-f)-X(g-x) = 0, 

i. e. the Weierstrass condition is satisfied. 

Insert into consideration a scalar function U — U(x, X,t) doubly continuously differen- 
tiable with respect to its arguments for the Hamilton system (1.9), (1.10), which we name 
the controlling function. That the function U(x, X,t) is supposedly not given beforehand, 
but found from quite definite target conditions. 

The controlled mapping from the old variables x(t) G R n , X(t) G R n to the new ones 
y(t) G R n , p(t) G R n is realized using the relation 

x + U x = y, X -U x = n, (1.16) 

where U\ = dU/dX, U x = dU/dx. Note that the record x = q, X = p, y = Q, p, = P is 
usual, further we however follow the established notation. In addition, it from expressions 
(1.16) follows that in case of identical transformations x = y, X = p, we have U\ = 0, U x = 
0, and U = may be taken as the controlling function U — U(x, X,t). 

Give attention to the fact that the mapping (x, X) — > (y, p) given by the equalities 
(1.16) is the generalization of the known Hilbert-Courant mapping [5] for the paramet- 
ric generative function T(a,b,t) in a canonical system (a and b are here the vector pa- 
rameters) at the change from the variables to q — q(a,b,i), p = p(a,b,t) to those of 
Q = Q( a , b,t), P = P(a, b, t), by the rule 

q = a- l -Y h , p = b+ l -Y a , Q = a + l - T b , P = b-^T a , (1.17) 

where T a = dF /da, = dT/db. The removing parameters a and b give in these formulas 
clearly give rise to the relations of form (1.16). 

An important advantage of mapping (1.16) over the parametric form of mapping (1.17) 
consists in the lack of any intermediate parameters and in the subsequent necessity for 
their determination. 

The Jacobian of mappings (1.16) is assumed to be non-zero: 

det (—) = D ^ yi ' — — ^ 
\dxj D(xi, x n ) 

det UJ = d(m, .... a.) * °' (L18) 

In relations (1.18) there are determinants of the matrices E ± d 2 U / dxdX respectively, 
where E is a single- valued matrix of n dimensions. 



Instead of relation (1.16), a controlled sympletic mapping may be taken that resembles 
by its form the Hilbert-Courant mapping 

x + ±IU x = y, \- l -IU x = ^ I= (^ E o)' (L19) 

where E is the sympletic matrix. It is easy to see that in the presence of (1.19), the 
Jacobians are the same and equal to 

det (l) =det (l) = 1 + I det (lsx)^ 

We however confine ourselves in our further arguments to mappings (1.16). 



2 Essential traits of controlled mapping 

The leading question that interests us at the given stage is what qualitative changes are 
introduced by the controlled mapping (1.16) into the description of canonical system (1.9), 
(1.10) 

x = H x , A = -H x 

with the Hamiltonian dependent on the canonical arguments x, A, t. 

We have known that the invariant Hilbert integral (1.14) calculated over the extremal 
transforms into the action functional S that takes a minimum value, see principle of the 
least action in the Hamilton form 



S — [ L(x,x, X,t) dt — [ Ldt, 
J to J(Co) 



independent of the integral path, but dependent only on the ends of this path. We have 
over the extremal 

S= [ Ldt= [ (^x-H]dt= [ Xdx-Hdt. (2.1) 

J (Co) J (Co) V OX J J (Co) 

The integral (2.1) is the above treated Hilbert integral. The differential expression Xdx — 
H dt being its member has the form of an integral Poincare-Kartan invariant. 

We assume that in the expanded phase space R 2n+1 with coordinates x G R n , A G 
R n , t G R and the Hamilton function H = H(x,X,t) constructed for the system (1.9), 
(1.10), the controlled transformation (1.16) is diffeomorphic, where the diffeomorphism is 
considered as a mutually single-valued and mutually differentiable mapping. 

Let us construct 1-form u 1 = Xdx-Hdt and refer to the following definition [10, 
p. 859]. 

Definition 2 Diffeomorphism retaining the external differential form u 2 = Yh=\ dXiAdxi 
is referred to as a canonical mapping. 



As the phase trajectories (phase flow) (x,X) of system (1.9), (1.10), representing [1] 
the rotor lines for the form Xdx — H dt such that by the Stokes lemma {§ oo 1 = § l2 u 1 ,) 
the statement (the theorem on the integral Poincare-Kartan invariant) follows that 



j> Xdx — H dt = j> Xdx — H dt, 

J 71 J 72 



where 71,72 are the closed curves enveloping one and the same tube of phase (integral) 
trajectories of Hamilton system (1.9), (1.10). 

In these relations the form A dx is referred to as a relative integral Poincare invariant. 
For the two-dimensional part of the tube of rotor a, the Stokes formula 

Xdx = j j dX A dx 



7 



is valid, the integral invariance then follows of the 2-form u 2 = dX A dx for the phase flow 
(x,A)[l]. 

The canonical transformations convert the Hamilton system into a Hamiltonian system 
as well. We in this connection turn our attention to a known theorem. 

Theorem 1 Let F : R 2n — > R 2n be the canonical transformation of phase space that 
converts the point with coordinates (x,X) into that with coordinates (y, /j,). In the new 
phase coordinates (y,/J>), canonical equations (1.9), (1.10) then have the canonical form 

dy dG d\i dG 

dt = fy? ~dl = ~l)y~ ( ' ' 

with the former Hamilton function 

G(y,fjL,t) = H(x,X,t). (2.3) 

This theorem is proved with the use of a canonical property of mapping, namely: 

<t> Xdx = <t> /j, dy 

J 7 J'y 

over any closed curve 7, from where 

f X dx — /j, dy = 

and, therefore, the integral f^l'xo) ^dx—jj,dy = Q, where J^xo) Xdx—^idy = Q depends 
on no path of integration, but an initial (x , A ) and final (27, Ai) point of path. Here we 
have 

Xdx — fidy = dQ. 

Q = Q(x,X,t) is here an arbitrary continuously differentiable function of its arguments. 
This equality in an expanded phase space R 2n+1 transforms into the equality 



Xdx — H dt = jidy — H dt + dQ. 



(2.4) 



The trajectories of canonical system (1.9), (1.10) are represented by rotor lines of the 
form fidy — G dt + dQ (see details in [1]). Comparing this form with form (2.4) in the 
right side, we therefore arrive at the conclusion that equality (2.3), G(y, fi, t) = H(x, A, t), 
is valid. If the Hamilton function does not change in the canonical transformation, i. e. 
G — H, such a transformation calls quite a canonical transformation. 

If we have the considered case of controlled mapping (1.16), y = y(x,X,t), \i = 
fi(x, A, t), using formula (2.4), where on the right there is G instead of H, we arrive at an 
important conclusion that in the canonical transformation of phase space dependent on 
time t, the canonical equations (1.9), (1.10) in the variables y,fJ>,t have a canonical form 
(2.2) with a new Hamilton function 

dO 

G(y,n,t) = H(x,\,t) + -^, (2.5) 
where, as before, we have (x , A ) at the fixed initial point 

*(xi,Ai) 
'(xo,Ao) 



Q(x\,\i,t)= / Xdx — \idy 

J fxn.An) 



So, for the Poincare-Kartan invariant to retain its form, we must require that the 
changed part of this integral invariant is a total differential: 

dQ dQ dQ 
Xdx — H dt = tidy — G dt + — — dx + — — dX + — — dt, 

dx dX dt 

from where we obtain the system of relations whose fulfillment must be required for 
providing the canonicity of the controlled transformation: 

G — H — Q t , (X - Q x ) dx = fidy + Q x dX. 

To use the results of Theorem 1, it needs to demonstrate what conditions are required 
for controlled transformation (1.16) to be canonical. 

Theorem 2 Let conditions (1.18) be fulfilled and, in addition, the controlling function 
U = U(x,X,t) and the variable X(t) G R n in transformation (1-16) satisfy the differential 
equality 

(U x -X)dU x = U x dX, (2.6) 

where dU\ = U\ x dx + U\\ dX + U\ t dt. The transformation (1-16) of variables x, X — > y, /j, 
will then be canonical. This transformation changes Hamilton system (1.9), (1.10) with 
the Hamiltonian H(x, A, t) into the Hamilton system (2.2) with the Hamiltonian G(y, /i, t) 
by rule (2.5), where Q(x, X,t) = U (x, A, t). 

Proof. The canonicity criterion similar to canonicity criterion given in [2, 3] actually 
follows from (2.4), (2.5): if in the transformation of variables x, X — > y,/i the differential 
form 

W = ndy-Xdx-{G-H)dt (2.7) 



is a total differential for a function V = V(x,X,t), i. e. W = dV, where dV = V x dx + 
V\d\ + V t dt, we have a canonical transformation. 

In the considered case, the arguments x,\,y,fi in the Hamiltonians H and G are 
related through expressions (1.16), and the difference G — H satisfies equality (2.5). We 
must in this way take the function — Q as a function V, following equality (2.4). Setting 
further V = — Q = — U, we obtain the following expression for the form W (2.6): 

W = d(-U). (2.8) 

Let us verify the fulfillment of this condition. Substitute relations (1.16) and (2.5) 
into the form W (2.7). We then get 

W = (A - U x ) d (x + U x ) -\dx-U t dt = (\- U x ) dU x -U x dx- U t dt. 

With regard for the theorem condition (2.6), we obtain for expression (2.8) 

W = - U x d\-U x dx-U t dt = - dU. 

It has thereby been argued that the controlled mapping (1.16) is a canonical mapping 
that establishes the theorem itself. 

Some arguments of a general nature, which are immediate from Theorem 2 that deals 
with the canonicity of the controlled mapping, can be proved. Note that the scalar 
condition (2.6) may be written in a form all the more compact 

U x dU x = d(\U x ). 

The condition (2.6) is in any case the equal of the canonicity criterion on the existence 
of some function V(x, X,t) such as W = dV, where the form W is described by equality 
(2.7). 

How could we take advantage of canonicity condition (2.6)? In principle, criterion 
(2.6) can be seen as the criterion that is satisfied by the controlling function U(x, \,t) (or 
rather its vector differential components U x , U x ) and vector-function of variables A(£). 

However, to limit the choice of the function U, on the assumption of condition (2.6), 
where the vector A is given by equation (1.9), is a rather strict requirement. Because it 
first of all groups together 2n unknown vector-functions U x and U x . The condition (2.6) 
would therefore be appropriate for the choice of a not uniquely defined vector A rather 
than in the choice of U. 

Now demonstrate how this can be done. Let the Jacobi matrix as a preliminary 
be denoted as f x = A, A = A(x,t). Then it is known [11] that the solution of vector 
differential equation (1.9) X(t ) = A may be written as 

\(t) = B(x,t)\ , 

where B(x,t) is the matrix obeying the matrix differential equation 

B(x, t) = - A(x, t) B(x, t), B(x , t ) = E, 

where xq = x(to) E is the unit n-matrix. The formula (2.9) for the mapping Ao — > A 
determines the differential homomorphism of class C 1 . 



Substitute expressions (1.9) and (2.9) into canonicity criterion (2.6): 



(X-U x )dUx = U x AXdt, 
from where we arrive at the scalar equation 

(U X -A*U X )\ = U X U X , 

or 

(U X -A*U X )BX = U X U X: 

where * over indicates the transposition operation and all vectors are multiplied as scalars. 
Let us also denote the vector C = B* (U x — A*U X ). Then the latter equation is 

CA = U X U X , 

where CXo = (C, Ao) = Yh=\ CjAoi may be considered as an equation with respect to one, 
for example, k-th initial condition Xot : 

_ U X U X — Yh-i C{\oi 

A 0k — ~ \t=to: 



where at the top of the formula there is i — 1, n, i ^ k, and to is the fixed initial moment 
of time. 

One more point needed to be made that the controlling function U, by the before made 
premises and especially the assumption that U = Q, U(xi, Xi,t) = f^Xo) Xdx — fidy, is 
of the nature of potential (energetic) function. In this connection the function U may be 
named the controlling potential function, or controlling potential. 

It is significant that the canonical condition (2.6) is realized with a feedback. If the 
initial moment of time to to be really considered as a fixed one, the expression for Ao& is 
determined, having regard to the solutions of the corresponding equations of dynamics 
(1.9), (1.10), i. e. to the action-functional extremals, through the current values of phase 
variables. We arrive by that at the problem of canonical synthesis as far as here we 
construct the law for forming the initial value Aofc, which provides the fulfillment of the 
canonicity condition for a controlled mapping and covers all manifold of other initial data. 

Let us show that the above obtained condition of canonicity and the Hamilton- Jacobi 
equation for the controlling function are closely related. 

Theorem 3 Let canonicity criterion (2.6) be fulfilled for mapping (1.16). Then, in 
order that mapping (1.16) to be a solution of the Cauchy problem for Hamilton equations 
(2.2) 

dG dG 
y=-QjJi ^ = ~~dy' yo = y(t ), /i = /i(t ), 

the controlling function U must satisfy the Hamilton- Jacobi equation with the Hamilto- 
nian G : 

_ =G ( I+ _ A __ ( ), (2 . 10) 

where U = U(x,X,t), U(x, A, to) = 0. 



Proof. To prove the theorem, we need equality (2.5) given at Q = U. If canonicity 
criterion (2.6) is fulfilled, we have 



dU_ 
~dt 



+ H(x,\,t)=G(y, l i,t). 



(2.11) 



It remains to be noted that in the Cauchy problem for the system of canonical equa- 
tions (2.2), the variables x and A became initial points of the canonical system 



for the trajectory y = y(t), /i = (lit) with the Hamiltonian G(y,fx,t). As this takes 
place, the motion equations are shaped into the simplest form that appropriates to a zero 
Hamiltonian H : H(x, X,t) = 0. Therefore, to determine the controlling function U from 
equation (2.11), we shall obtain equation (2.10). This completes the proof. 

It is ready to discover that for the Cauchy problem of system (1.9), (1.10) at x = 
x(t ), A = A(t )- Theorem 3 can be reformulated to terms of an old Hamiltonian H(x, A, t) 
for Cauchy's problem of system (1.9), (1.10). We in this case have y = x(t ), \i = A(t ), 
and G(y,/i,t) = in equation (2.11). 

Theorem 4 For the Hamilton system (1.9), (1.10), the Cauchy problem solutions may be 
presented by canonical mappings (1.16) when the function U satisfies the Hamilton- Jacobi 
equation with the Hamiltonian H having the form 



where, as before, U = U(x,X,t), U(x,\,t ) = 0, and the variables x and A are solutions 
of equations (1.9), (1.10) respectively. 



(1) We can see from the structure of the proofs that Theorems 3 and 4 are indeed 
invertible, i. e. represent the necessary and sufficient conditions for the solutions 
of corresponding Hamilton-Jacobi controlling function equations (2.10), (2.11) to 
exist. 

(2) Draw our attention to the fact that the functions U — U(x, A, t) are common in 
their notation in equations (2.10) and (2.11), but different in values. 

(3) The controlling function U should not be confused with the generating function 
that is actively used in the classical formalism of the theory of Hamilton equations, 
canonical transformations, and integration of differential equations by the Hamilton- 
Jacobi method. 

These functions are different in meaning of their formation and further application. 
Their main distinction consists in the dependence of the controlling function on the old 
variables x, A only, including time t. Recall that the generating function U is an arbitrary 
function of mixed (old and new) variables. 



x = y(t ), 



A = //(t ) 



dU_ 
~dt 



+ H(x,X,t) = 0, 



Remarks 



Of course, the function U can be given the nature of generating function U . With this 
aim in view, let us write mappings (1.16) in a general form, 

y = ip(x,X,t), ^ = ip(x,X,t), (2.12) 

where f = x + U\, if) = A — U x , assuming that at d 2 ip/dxdX ^ 0, d 2 ip/dxdX ^ they 
are solvable in the old variables 

x = x(y,fi,t), \ = \(y,n,t). (2.13) 

The substitution of one of formulas (2.13) into the function U for x or A (there can 
be nothing but four variants) leads to the appearance of a generating function U . The 
advantage in the use of controlled mapping (1.16) with the controlling function U over 
the mapping with the generating function U is obvious: the use of U does not provide for 
the resolution (reversibility) of equations (2.12) in x and A and for the conversion to the 
explicit dependences (2.13). 

It was specified in [1] that the generating function formalism seems to be 'depressive 
in its non-invariance and essentially uses phase-state coordinate structure'. Taking into 
account this note, let us study the question on invariant properties at controlled mapping 
(1.16): (x,\,t)^(y,fi,t). 

Consider, for example, the converted form of the integral Poincare-Kartan invariant 
in the canonical conversion (x, A) — > (y, /i) under the generating function U (x, //, t) 

\dx- Hdt = fidy-Gdt + dU, (2.14) 

where dU = U x dx + C/ M d/i + U t dt. It follows in particular that 

G = U t + H, X = U X , y = U IM . 

The canonical transformation (x, A) — > (y, fi) obtained with the generating function 
U (x, /j,, t) is admissible if the condition 

det ° 2 - ^ 
dxdfi 

has been fulfilled in relation (2.14). It is obvious that this condition depends on the choice 
of new canonical variables. 

On the contrary, for the controlling function U(x, X,t) with the converted form 

Xdx - H dt = fidy - G dt + dU, 

the existence condition for the canonical transformation directly depends on the condition 
det (E ± d 2 U/dxd\) ^ (1.18) applied to the controlled mappings (1.16). This condi- 
tion is invariant with respect to the new canonical variables and may be affected by the 
canonical change of variables. 

Consider further the Lagrangian K derived from the Lagrangian L(x, x, A, t) by means 
of controlled change of variables (1.16) if we set in addition that relations (1.16) are 
resolvable over the old variables x and A by formulas (2.13): 



L(x,x,X,t) = L[x(y,/i,t), x{y,n,t), X(y,/i,t), t] 



= K(y,y,fi,fx,t), (2.15) 
where x(y, /i,t) = x y y + x^jj, + x t . Here we have 

dx dx dx 

y = ^' ""^d^' Xt= dt 

in standard notation. 



Definition 3 The Lagrangian K (2.15) obtained after the change of variables in the 
Lagrangian L, will be called an induced Lagrangian. The corresponding Euler equation for 
K will be called the induced Euler equation and its solution the induced extremals. 

The following theorem on the invariance of mappings (1.16), (2.13) is valid, which 
generalizes a known statement of calculus of variation [9]. 

Theorem 5 Let all suppositions about the existence of mappings (1.16), (2.13) be ful- 
filled. Then the extremals for the variational problem with the Lagrangian L are the 
induced extremals for the variational problem with the induced Lagrangian K. 

Proof. Let us use the notation v = K y , i. e. 

OK dL dx . _ /n „ „. 

v = IT = 7T7T = Afi > 2 - 16 
ay ox dy 

as far as x y — x y . We find now K y : 

dK dLdx dL dx dL d\ (2 17) 

dy dx dy x y d\ dy 

= L X Q + \E + (x- /)A, 



where the notation 



dx d\ 
dy' dy 



has been used. 

Let us find the value of matrix S from expression (2.17): 

„ dx d . . d 2 x d 2 x d 2 x 

- = tt- = Tr{x y y + XpV + Xt) = TT^T V + a a A* + 



dy dy dy dy dy d/i dy dt 

On the other hand, the value of matrix Cl, where Q = x y = £l(y, fx, t), may be written as 

• dQ dQ dfl d 2 x d 2 x d 2 x 
" = -zrV ~i k A 4 H ttt = tt^t v + ^r^r A* 



dy dfx dt dydy dfidy dtdy 
It from here follows that we have the equality E — Cl. 



Then we construct an induced Euler equation with the induced Lagrangian K using 
equations (2.16), (2.17): 



= = S («l)-Wl-«l + (*-/)A 

= (A - L x ) Si + (x - f) A. 

It thus follows that the equations A = L x , x — f, and consequently v — K y , take place 
on the extremals of the Euler equations generated by the Lagrangian L. That establishes 
the theorem. 

Note in passing that we would arrive at the conclusion which has been formulated 
in Theorem 5 that mappings (1.16), (2.13) are invariant if we use the invariance in the 
Poincare-Kartan form when the choice of variables to be canonical. 

The representation of the controlled mapping is, of course, not be limited by formulas 
(1.16). They were entered for definiteness sake and for the demonstration of the relation 
with Hilbert-Courant mapping (1.17). The same treatment may be applied very well to 
the controlled mapping of the form 

x±U x = y, X±U x = pi, (2.18) 

or 

x±U x = y, X±U x =fi (2.19) 

may be as readily entered, the signs on the left parts of equations (2.18), (2.19) can be 
either identical or opposite. Still more exotic combinations seem to be possible when 
constructing new variables y, \i. 

The main requirement for the new variables, nevertheless, as well as for the old ones, 
is their canonicity. Take, for example, a controlled mapping of the form 

x + U x =y, X — U\ — fx. (2.20) 

Satisfy that the transformations (2.20) lead in specific situations to the canonical 
variables y, //. In this theorem analogous to Theorem 2 may be proved. 

Theorem 6 Let conditions (1.18) be fulfilled with the controlling function U = U(x,X,t) 
such that its partial derivatives U X ,U\ satisfy equations (2.20). If the equality 

(A - U x ) dU x = (U x - U x ) dx - U x dX (2.21) 

is then satisfied, the transformation of the variables x, X — > y, /i will be canonical. System 
(2.2) has here the Hamiltonian G(y,/i,t), (2.5), where Q(x,X,t) = U(x,X,t). 

Proof. Theorem 6 schematically appears in the following form. We formulate differential 
form (2.7), where W = dV, V = V(x, A, t) is any function. Further it needs to check that 
if the canonicity criterion (2.21) is fulfilled, the function V = —U, where U = U(x,X,t) 
is the controlling function, it is a desired one. We have by mapping (2.20) for form (2.7) 



W = (A - U x )d(x + U X ) - Xdx - U t dt 



= XdU x - U x dx - U x dU x - U t dt 



= (A - U x ) dU x -U x dx- U t dt 

= - U x dx -U x d\-U t dt = d (- U). 

The latter record is, obviously, completes the proof of the theorem. 

In addition, the canonicity criterion (2.21) can be provided by the choice of a certain 
k-th initial condition A fc. In this connection to solve equation (1.9), we take equality (2.9). 
Considering initial condition (1.1), criterion (2.21) can be written in the following form: 

(A - Ux) dU x = [(Ux - U x ) f + U x A\]dt, 

hence, we get the equation 

(A - Ux) U x = (Ux -U x )f + U X A\, 

or 

(U x - A*U X ) B\ = U X U X + (Ux - U x ) f. 

An expression for A fc remains to be written. Resolving the latter scalar equation in 
A fc, we obtain the formula 

, UxU x + (Ux-U x )f-Z? =l DjXoi 

A 0k — p: , 

Dk t=t 

where % ^ k, k — T~n, D = B* (U x - A*U X ). 

Example 1 The controlled mapping (2.2) was chosen not accidentally. It is convenient 
with its help to make rotations of the phase coordinate system (x, A). 

For example, the rotation of axes x and A through the right angle in an anticlockwise 
direction that is executed by transformations (2.20) must be consistent with the equalities 
(y — A, (j, — — x) : x + U x = X, X — Ux = — x, from where the expressions for U x and Ux 
follow: 

U x = X — x, Ux = X + x. 

The controlling function 

A^ x^ 

U = U(x,X,t) = — + Xx + u(t), A 2 = (A,A), x 2 = (x,x) 

obviously provides an appropriate rotation of axes x and A for every continuously differ- 
entiable function u(t). 



3 Energy transformations and controlling fields 



By the controlling field is meant a scalar functional field in R that is given by the values of 
function U(x, X,t). We previously mentioned that the controlling function U — U(x, X,t) 
was of energy nature from the assumptions made on the canonicity of a controlled map- 
ping. Turn to this point. 

Take once more expression (1.2) for the action functional S. The action functional, 
where t — t\ in integral (1.2), t being the running time, and the initial point (xo,Xo,to) 
fixed, is said to be an action function 



S(x,X,t)= f L(x,x,X,t)dt. (3.1) 



Try to discover a relation between the action function S (3.1) and the controlling 
function U. If such a relation will be estimated, we shall be able to form an algorithm for 
determining the controlling function. 

Denote through 7 the extremal joining its initial point (x , t ) with its terminal at the 
running point (x,t). For the action function S(x, X,t) (3.1), we thus have as well 



S(x,X,t)= / L(x,xX,t)dt. 

J-y 



According to [1], we may demonstrate that if the initial point (x , t ) is fixed, the differ- 
ential of the action function S(x, X,t) has the form of integral Poincare-Kartan invariant 
(cf. with expression (2.1)): 

dS = Xdx- Hdt, (3.2) 

where 

, n dS , OS „ OS , x dL 
db = — dx+ — dX + — dt, A = — , 
ox OX at ox 

and the quantity H = Xx — L is determined in view of a finale velocity x of the trajectory 
7- 

The proof of equality (3.2) is based on lifting the extremal 7 from the space (x,t) 
to the expanded phase space (x,X,t), where A = dL/dx. In this case, the extremal is 
substituted by a phase trajectory that is among the variety of rotor lines having the form 
Xdx — H dt. 

By the fact that the relations 

dS „. , x fdS \ , dS „ 

— = -H(x,X,t), ^—-X)dx + —dX = (3.3) 

follow from equality (3.2), we conclude that the action function S = S(x,X,t) satisfies 
the Hamilton- Jacobi equation 



dt \ dx 
The form (1.2) can indeed be written as 



dS = L dt, 



or S = L, where L = L(x, x, A, t) — X (x — /). 
We have then 

S x x + S x \ + S t = L, 



from where the equalities 

S x — L± = A, 



0, S t = L-L x x = -H = -Xf 



follow (compare with equalities (1.13)). With an additional condition S\ — 0, the first 
equation in system (3.3) is therefore converted to form (3.4). 

Having regard to the written relations, the Hamilton- Jacobi equation (3.4) of the 
scalar action function S(x, A, t) can then be presented in the form 



as_ _dS_ 

dt dx 



f(x,t). 



The system of 2n equations 



— - A — - 

dx dX 



has in this case to be thought of as a system of equations in In unknowns x(t) and X(t). 

If the Cauchy problem with the initial condition S(x, A, to) — Sq{x, A) is formulated for 
equation (3.4), its solution is then reduced to the solution of canonical Hamilton equations 
(1.9), (1.10): 

x - dH X 8H 
dX ' dx 

with the initial conditions 



xito) = x Q , X(t ) = X 



dSo 
dx 



(3.5) 



X=XQ 



The solution of this problem on the space (x, t) comprises, as we know, the extremal 
x = x(t) for the given variational principle A / Ldt = 0. This extremal issues out of the 
initial point x and is called the characteristic of problem (3.4). We integrate equality 
(3.2) along the characteristic that joins the points A = (x ,t ) and A = (x,t). And for 
the action function S with the initial condition S , we find the expression 



S(A,X) = S (A ,X ) 



L(x, x, A, t) dt, 



(3.6) 



which gives the solution of the problem (3.4). Note in addition that the initial condition 
Ao must be correlated to the initial condition Aofc (3.5). 

Let us briefly consider the most important autonomous case. Let time t now not 
come explicitly in the expression of Hamilton function H, i. e. OH /dt = 0. The equality 
dH/dt = dH/dt was above obtained in view of he Hamilton equations (1.9), (1.10). We 
therefore have dH/dt = in the given variant, from where it follows that H = h = const 
is the first integral (Jacobi integral) of equations (1.9), (1.10). 



By the Legendre transformation, H = (dL/dx)x — L = Xf and / = f(x). The 
Legendre function L is then independent of time too; L = L(x,x,X) = X (x — f) and 
dH/dt = - dL/dt = 0. 

Let the surface H(x,X) = h be projected from the expanded phase space (x,X,t) to 
the space (x, A). In this case, the time t actually does not vary ((dt = 0)), the total 
variation A changes into the isochronous variation 5, and the expression (3.2) takes a 
shortened form as compared to the relative integral Poincare invariant: 

dS = Xdx. (3.7) 

Hence, the trajectories of Hamilton system (1.9), (1.10) are extremals for the variational 
principle corresponding to form (3.7). 

The phase trajectories of canonical equations (1.9), (1.10) lying on the surface H(x, A) = 
h are in this way extremals of the integral 



= f Xdx, (3.8) 



and these extremals joint the points xo and x\. The formulated principle forms the 
contents of the Maupertuis-Lagrange least (stationary) action principle validated by La- 
grange [12] with regard to the Lagrange action J, namely, 



j-t'2 

J = / 2Tdt, 



where T = mv 2 /2 is the kinetic energy of a material point of mass m and velocity 
v = ds/dt, namely 

J = mv dt = mv — dt — mv ds. 
h x h x dt Jsi 

It is clear that in the agreed notation the Lagrange action J will be equal to the 
shortened action integral S (4.8) if the impulse mv = X and distance s = x are taken. 

Let the controlled mapping (1.16) be in the autonomous case given by the controlling 
function U = U(x,X). If this mapping is quite a canonical one, then by Theorem 1, 
the conversion from Hamilton equations (1.9), (1.10) to those of form (2.2) takes place 
according to equality (2.3) for the Hamiltonians G = H, where G = G(y, /i), H — H(x, A). 

In so doing, the relation Xdx — jidy = dQ is valid (see expression (2.4), where Q = 
Q(x, A)). We have by that 

Xdx = n, dy + dx + -Q- dX, 
ox oX 

we get from where 

(A — Q x ) dx = \idy + Q\ dX. 

Select as a function Q : Q = U, U = U (x , A) . Then for the canonicity of transformations 
(1.16), the equality 



(A - U x ) dx = (A - U x ) d(x + U x ) + U x dX 



needs to be required, or 

d(XU x ) = U x dU x , 

i. e. we again arrive at the canonicity criterion (2.6). 

We shall return to basic relation (2.4), where Q(x,X,t) = U(x,X,t) and the func- 
tion U(x,X,t) specifies controlled mappings (1.16). These mappings are canonical when 
criterion (2.6) is fulfilled, i. e. at a certain choice of the initial vector Ao- 

Further, let us use the concepts of Lagrangian, Hamiltonian, action function, canon- 
ical variables, and write energy the equality (2.4) being main in the controlled mapping 
method in the form 

Ldt = Kdt + dU, (3.9) 

where L dt = Xdx — H dt, L = L(x,xX,t) is an initial Lagrange function, x, X are ini- 
tial canonical variables, H = H(x,X,t) is the initial Hamilton function, K dt = fidy — 
Gdt, K = K(y,yfi,t) is a new Lagrange function: 

K = fi(y-g), V = Ky, G = fig. (3.10) 

G = G(y, fi, t) = fig is here a new Hamilton function, where g = g(y, t), which is connected 
with the function if by a Legendre transformation and y, fi are new canonical variables 
in the new (2n + l)-dimensional expanded phase space (y,fi,t). 

The function K (y, y, fi, t) is corresponded by a new action function R = R(y, fi, t): 

R(y,V,t)= I K(y,y,fi,t)dt= f K(y,y,fi,t)dt, (3.11) 

where S is the extremal of new variational principle A / K dt = 0, 5 being an induced 
extremal. The new canonical equations with a new Hamiltonian G have the form of 
equations (2.2) 

• dG dG 

where g = g(y,t) is a certain function continuously differentiable in its arguments. 

It is important to note that the choice of the new Hamiltonian G = G(y, fi, t) and of the 
right parts in canonical equations (3.12), respectively (for example for a more convenient 
integration or reduction of the systems), in a more simple form (its normalization) is a 
main purpose condition of the formation of the controlling function U — U(x, X,t). 

Using the formalism of action functions, we write relation (3.9) that determines the 
controlling potential U (x, X, t) in the form 

dS = dR + dU, S-R = U, (3.13) 

from where we derive the expression for U(x, X,t): 

U(x,X,t) = J Xdx- Hdt- (jf fidy-Gdtj (3.14) 
= / Xdx- Hdt- / fidy-Gdt) = (Xx - fiy + G - H) dt. 



In formula (3.14) by the canonicity of controlled transformations, the written line 
integrals do not lye on the path of integration but depends only on the finale values 
(x,X), (y,/i) at the fixed initial values (xq, Ao), (yo,A*o)- 

In addition, in (3.14) the new canonical variables y,/i are connected with the old 
canonical variables x, A by controlled mappings (1.16). Note also that for the new action 
function R (3.11) (from equation (3.13)) there takes place a Hamilton- Jacobi equation of 
the form 

Rt = — G, Ry — /i, R^ = 0, 

where G = G(y, R y , i). 

Example 2 Consider an autonomous dynamic system that is square integrable and de- 
scribed by the equation 

x = f(x), x £ R n 

with the given initial condition x = x(t ) and suitable requirements to the vector-function 
f(x). It is now necessary, using controlled transformations of the form (1.16), to reduce 
it to the given system presented by the equation 

y = g(y), y e R n , 

where the vector-function g(y) is determined from the given Hamilton function G(y,fi). 
Use all previous notation and constructions to solve this problem. We have 

L(x, x, A) = A (x — /), H(x, A) = Xf = h, 

where h > is a constant equal to the value of total system energy. Canonical equations 
in the variables x, A have the form 

i = ^ = /(*), X = -^ = -f x (x)X. 

It is considered that the initial data x Q , A are given and A is coordinated with the 
canonicity criterion (2.6). 

Let the new Hamiltonian G meet the requirement 

G(y,fi) = a/i, 

where a G R n is a given constant vector. Due to the autonomy, H = G = h. The new 
canonical equations (3.12) may be written in the form 

. dG dG n 

y = j— = a, n = - — = 0, 
oil oy 

from where it follows that y(t) = at + b, ji(t) = c. Here a, b, c are the vectors of known 
components. Thus, we have g(y) = a, G(y,/i) = ac = h. 

It is required, starting from these purpose conditions, to determine the controlling 
function U — U(x, A) such that it is reduced to these new Hamilton equations. Let 
there be controlling mappings (1.16) and by the choice of Ao) canonicity criterion (2.6) 
be provided. 



Here with regard to (3.9)-(3.14), we have for the situation of finding the controlling 
function U (x, A) : 

\dx = [idy + dU, U = S — R, 
from where, since H = G = h, we get a chain of equivalent relations 

U(x,X) = J Ldt- jf Kdt 
= X[x — f(x)] dt — c (y — a) dt 

J 8 

r(x,X) r(y,c) Ax,\) 

= / Xdx — cdy — Xdx — c(y — y ), 

J(xo,Xo) J(yo,c) J(xo,Xo) 

where y = y(t), yo = y(to), and by the canonicity of transformations, all written line 
integrals do not depend on the path of integration, but depend only on the finale values 
at the fixed initial data. 

Substituting the expression for y (1.16) into the latter equality, we obtain an equation 
for the required function U(x, A): 

r{x,\) 

U(x,X)— / A dx — c (x + U\) + cy . 

•A^o.Ao) 

This equation is a first-order partial linear differential equation with respect to the 
unknown function U(x, A). This has the form 

U(x,X) + c^^ = F(x,X), 

where the notation is 

dU(x,X) " dU(x,X) XN f(.*M , , 

F(x, A) being a known function of variables x and A. 

An attempt to get an analytical solution of equation (3.15) meets certain difficulties 
[13, 14] caused by 2n independent variables xi, x n , X±, A n and n partial derivatives 
U\ 1 , U\ n . To obtain the relation between x±, x n , Ai, X n and U, we must find a 
solution for the corresponding system of ordinary differential equations (characteristic 
equations of a specific form) that satisfy the given initial conditions. If we disregard the 
theoretical questions concerning the solution of equation (3.15), we can point to the fact 
that in the given case this equation can be considered as a final algorithmic equation for 
determining the controlling potential that gives a mapping in the initial phase space (x, A) 
with subsequent numerical resolution. 

Let us apply the obtained results of the method of dynamical controlled-mapping 
method to the most important class of mechanical systems, Hamiltonian systems, which 
in standard notation of generalized coordinates q(t) G R n and generalized momenta p(t) G 
R n , have the form 

dH dH 



3.15 



where H = H(q,p,t) = pq — L(q,q,t) is the Hamilton function, L = L(q,q,i) is the 
Lagrange function, L — T — II, T is the kinetic energy of the system, II = Tl(q, t) is 
its potential energy, T = T2 + T\ + To, T2 = (1/2) X) dijqiqj is the quadratic form of 
generalized velocities, A = (a^) is the kinetic energy matrix, Ti = ^ is a linear form 
of generalized velocities, and T = T (q,t) is its null- form. Thus, we have 

j n n 

£ = ~ XI a *i + XI & * * + T o - n, 

2 ij=i i=i 

where the coefficients a^^bi depend on the generalized coordinates and time. Hence, the 
generalized momenta are expressed by 

Pi = jp- = 2^ aij^j + h- 
CQi j= i 

Since the kinetic energy matrix A is non-degenerated, solving then the latter linear 
equation with respect to q, we get 

n 

3=1 

where A' 1 = (a^). If to substitute this expression into the expression for H, we shall 
find 

n n 

H = X Pi9i - L = J2 a ij (Pi ~ bi) Pj 

i=l i,j=l 

1 n 

- 2 51 a ij a ik a ji (Pk ~ b k )(pi - h) 

i,j,k,l=l 
n 

- J2 b i aij (Pj ~ bj) - T + n 

j n n 
jj'=l i,j=l 

In case when the Lagrange function does not depend on time (dL/dt = 0), we have 
the first integral (generalized energy integral, Painleve-Jacobi integral) as follows: 

A . dL 

l^ViTrr- L = const, 
ti % 

or T 2 — T + n = const. If the system is conservative, T\ = T = and the generalized 
energy potential is the same as the total energy T2 + n = const (the total energy 
conservation law in the conservative systems). 

The Hamilton equations (3.16) are equivalent to the Lagrange equations 



where dL/d<ji = Pi are the generalized momenta, dL/dq,i are the generalized forces, the 
trajectories of motion for mechanical systems (3.16), (3.17) agreeing with the extremals 
of the action functional J t 1 L(q, q, t) dt (principle of least action in the Hamilton form). 

Further, we must rewrite the dynamic relations of mechanical systems in terms of 
customary standard notation x(t) and A(t). Let us take for this aim x(t) = (q(t), q(t)) or 
x(t) = (q(t), q(t)) and attach the meaning of vector Hamilton multiplier to \(t). We then 
resolve equation (3.17) with respect to a higher derivative and introduce the vector of new 
derivatives x(t). Now we arrive at a normal formulation of the mechanical system in the 
form of equation of motion (1.1). The neat process of controlled canonical transformation 
of system is completely corresponded to the above described scheme. 

It should be particularly emphasized that the controlled mapping method is based on 
the representation the system of equations of motion and Euler equations as a normal 
system of differential Cauchy equations (1.9), (1.10) respectively, giving the Hamilton 
function if in a special form (1.5), H = \f. Such a formulation sometimes leads to 
problems. 

Example 3 In the problem dealing with the determination of Lagrange multipliers for a 
ballistic flight (it is a section of path, where an engine is excluded) in a central gravitation 
field [7], it is necessary to present the equations of motion and Euler equations (Euler- 
Hamilton ones) in a normal Cauchy form and rewrite the Hamiltonian in the form of 
(1.5). 

The trajectory of motion, as we know from theoretical mechanics, is here represented 
by a conic section that in polar coordinates r, <p satisfies the equation 

P 

T = 

1 + ecos (</? — u) ' 

where we have the notation: p is a focal parameter ((p = r at <p = n/2 + u>), e is the 
eccentricity, u is a polar angle periphery center, i. e. the trajectory point closest to the 
focus of conic section, the angle a = (p — u being said to be a true anomaly. The change 
of polar angle is determined by the area integral 

r 2 <P = o- y/p, 

where a = \fyM, 7 is the gravitation constant, M is the mass of a central body. 

In order to compile the Euler equations for the multipliers A, we must first write the 
equations of motion in the normal form. Let us write them in polar coordinates. Use some 
known kinematic expressions for the projections of point acceleration on the generalized 
axes of curvilinear polar coordinate system 

.. . 2 I d(rV) 

w r = r - rip , w lf = — . 

r dt 

Since the point of mass m is under the action of only Newton force of attraction to 
the fixed center F = — mcr 2 r/r 3 , we obtain the equations of motion 

. 2 d(r 2 y) 



Note simultaneously that the second equation leads to an area integral, and the first one 
in passing to a new variable 1/r and new argument ip instead of t leads to the already 
known equation of conic sections. 

Of course, the equation that describes the motion of a point in a gravitation field in 
polar coordinate can also be found with the second-order Lagrange equations. For this 
aim, we must write the kinetic equation T and the generalized forces Q r and Q v : 

T = im(f 2 + rV 2 ), Q r = -m(^j 2 , Q v = 
and then construct the equations 

d or or „ 

dtW i ~W l = Qq - * = 1 ' 2 ' gi = r ' q2 = Lp 

being coincident with the above cited equations of motion. 

Construct now the equations of motion in a form solvable with respect to derivatives, 
i. e. found in their normal form. To this end, take the projections of the velocity of a 
point on the polar axes, namely v r = r, v v = rip. as velocity phase coordinates. As a 
result, we obtain the system of four equations of motion 

v l fcr\ 2 . v r v v v v 

v r = — ~ - , v v = S r = v r , f = —. 

r \r / r r 

The Hamilton function H (1.5) then takes the form 



and the system of Euler equations for the determination the multipliers X(t) will be written 
in the form 

-^1 = ^2 ~ — -^3; A2 = — 2Ai — + A2 A4 — , 

r r r r 

A3 = Ai — ^ — 2Ai —5- — A2 — 7T- + A4 — ^, A4 = 0. 

iy £ rjriQ tr*£i rjriZi 

The latter formula can be represented in the form A4 = const. The latter three form a 
system of linear equations in Ai, A 2 , A 3 with alterable coefficients. 



4 Canonical controlled mappings 

The controlled mappings (1.16) are constructed on the base of additive occurrence of 
controlling functions. The efficiency of their use may change from task to task. Point out 
in this connection that the transformation by formulas (1.16) can be of a more general 
nature being independent of the choice of initial data. 
Along with controlled mapping (1.16) 



y = x + U x ,n = X-U X ,U = U(x, A), 



let us consider a quasi-canonical transformation such that its vector variables x and A 
change by an infinitesimal value. Such a transformation is called the infinitesimal canon- 
ical transformation [1, 15]. Suppose 

y = x + rj(x, A)e,/i = A + q(x, A)e, (4.1) 

where rj(x, A), q(x, A) are the vector n-dimensional functions of x, A, and e is a small pa- 
rameter. If r) and q are given arbitrarily, the transformation (4.1) is clearly non-canonical. 

We therefore set forth the requirements imposed on the vector-functions f] and q that 
provide the canonicity of transformation (4.1), and find a formal analogy between the 
vector-functions U\,rj and U x ,q. 

Let us fall back on equality (2.4) of the form 

Xdx = fidy + dQ, 

where Q = Q(x, A) is an arbitrary function of x and A. Substitute (4.1) into its right side: 

A dx = (A + eq)(dx + e drj) + dQ. 
Dropping the terms that contain e 2 , we obtain 

e(q dx + A 77) = — dQ, 

or 

(q dx — 7] dX) + d(r)X) = — dQ. 
The latter equality may be written in the form 

q dx — r] dX = —d£l(x, A), 

where the notation 

£l(x, A) = -Q + Xr] 
e 

is entered for the scalar function Q(x, A). 

With regard to the fact that dQ = Q x dx + Q\ dX, we obtain from (4.2) 

q = —Q x , i] = n x 

due to the independence of the quantities dx and dX. Hence, canonical transformation 
(4.1), being infinitesimal, will be seen as 

y = x + fi A <M 11 — X- VL x q, (4.3) 

where Q = Q(x, A) is an arbitrary scalar function of variables x, A, which plays a great role 
(compared with (1.16)) of a controlling function with a small parameter or, by another 
formalism [16-18], of a small controlled perturbation (test signal). 

As for a general non-stationary case, all canonical transformations according to S. Lie, 
see [15], might be defined with one differential condition. Namely, the reversible trans- 
formation x, A — > y,fi (in this case the functional determinant d (y, fj,)/d (x, X) 7^ 0) is 
canonical if we have the differential form 



A dx = n dy + H dt + dQ, 



where H and Q are arbitrary functions of 4n + 1 variables x, X, y, /i, t, are satisfied iden- 
tically. 

This form should be given as another form that is well known if to subtract the quantity 
H dt from both parts of equation. Then we shall have 

Xdx- Hdt = fidy-Gdt + dQ, G = H - H . 

It should be noted that independent in these forms will be only 4n+ 1 variables x, X, y, /x, t 
from 2n+l because there exist 2n ratios between the variables of the form y = y(x, A, t), /i - 
/j(x,X,t). 

If we apply the canonical controlled transformation (4.3), where U\ = U\(x,X,q) = 
fl\(x, A) q, U x = U x (x, A, q) = Q x (x, A) q, to a scalar function ip(x, A), then it is possible 
to find its variation correct to the magnitudes of the order of q 2 smallness: 

ip(x i: X { ) = i/j(yi - U Xi , ^ + U x .) = ip(yi - Vt x .q, & + n x .q) 

" / dj) on | dip on 

~i\ dxj d\ d\i dxi 



or in the notation of Poisson bracket 

• dip dn dip dtt ' 

we obtain the formula 



~i\dxidXi dXidxi) 



ip(x u Xi) = ip(yi, m) - { ip, n } q = ip(y h Hi) - { ip, U }. 

Interest is provoked by the case of small canonical controlled transformation when we 
take the Hamilton function H(x, A) as a function Q and infinitesimal interval of time dt 
as an arbitrary small parameter. Suppose we have 

U(x, A, q) = Q(x, X)q = H(x, A) dt. 

For the infinitesimal canonical transformation (4.3), we may then write 

y = x + H\ dt, n = X — H x dt 

from where it follows by canonical Hamilton equations (1.9), (1.10) that 

dx = H\ dt, dX = — H x dt 

according to Hamilton canonical equations (1.9), (1.10) 

dx = if a dt, dX = —H x dt. 

These relations actually mean that in a time dt the system (x(t), X(t)) will pass to a new 
state (x(t + dt), X(t + dt)), and this change will be induced by the infinitesimal controlled 
canonical transformation such that the basic controlling function of the form 

d 

Q(x, A) = — U (x, X, e) 
oe 



is the Hamilton function H(x,X). We here assume that the initial function U linearly 
depends on the parameter q. 

In this way, it can be concluded, starting from this analysis, that the movement of 
this system is a sequence of infinitesimal controlled mappings, where the basic controlling 
function (at a small parameter q) must be the Hamilton function itself. Such a result only 
confirms an important conclusion [19] that the system can be led to the given specific 
state by means of sequential small probe signals (controlled perturbations). 

Turn to basic controlled transformations (1.16). We said in the preface to Section 4 
that it was possible to construct canonical controlled mappings without recourse to spe- 
cially selected initial data, as was done before. 

Theorem 7 Suppose the conditions (1.18) are fulfilled for transformation (1.16). Trans- 
formations (1.16) of the variables x, X — > y,/i are canonical if simultaneously : 

(i) to choice the vector U x from the equation 

U X = A*U X , (4.4) 

where A = A(x,t) = f x (x,t); 

(ii) to choice the vector U x using the orthogonality of vectors U x and U\. 

Proof. We shall, as before, start from a key equality (2.6) to which controlled mapping 
(1.16) leads. If the equality is fulfilled identically, the transformation (1.16) is canonical. 
Using initial Hamilton equation (1.9), we succeed in writing equality (2.6) in the form 

(U X -A*U X )X = U X U X . (4.5) 

The relation (4.5) will be fulfilled identically if we choose the vector U x such that equation 
(4.4) is valid, from where it automatically follows the requirement that the equality 

U X U X = (4.6) 

to be fulfilled. 

The vector differential equation (4.5) is identical in its structure to (1.9) A = — AX. 
The solution of equation (4.5) is given by the formula 

U x = D(x,t)U X0 , U X0 = U x \ t=tQ , (4.7) 

where D(x,t) is the n x n-matrix that satisfies the matrix differential equation 

D(x, t) = A*{x, t) D(x, t). D(x , t ) = E. 

The theorem is completed. 

In particular, it follows from (4.7) that the vector-function U x is a function of variables 
x and t. This means that the controlling function U (x, X, t) is linearly dependent on the 
A-vector elements: 

U{x, X,t) = / D(x,t) U X o dX + u(x,t), 

where u(x,t) is an arbitrary scalar function of x, and t that can be chosen, for example, 
from the orthogonality condition (4.6). 



Conclusions 



The idea to reduce a dynamic system to an expanded Hamilton form with simultaneous 
attachment of the formalism of classical calculus of variation was realized before as well 
(see, for example, [6-8] and [20, 21] for adaptive systems). 

It turned to be that such structures were able to work in the construction of controlled 
canonical mappings in Hamilton systems for the purpose of writing the equations of motion 
in a more comfortable form. 

For the purpose of constructing canonical coordinate mappings in Hamilton systems, 
the controlling-function method was developed in the given article, whose leading property 
was the dependence on the old canonical variables. 

This property is, in the end, allowed additional possibilities in forming an algorithm 
for a canonical change of variables, which render it redundant and less well defined as 
compared with the algorithm based on the generating-f unction method. 

Note also a broad range of new variational results, which, being of a general theoretic 
importance, were arranged as some theorems. These results may have practical use in 
the solution of various problems of canonicity and study of dynamic Hamilton- Jacobi 
relations. 

Notation 

R n is n-dimensional Euclidian space. 

Axi, 5xi are total and isochronous variations of the component Xi respectively. 
AS* is a total variation of action functional S. 
A* is the transpose with respect to the matrix A. 

S{C) = f^Qj L(x,x, X,t) dt is the integral S along the curve C, where x = x(t) is the 

representation of curve C,t e [to^i]- 

{■0, Q} is the Poisson bracket with the elements ip and fl. 

B x = dB/dx is the vector (matrix) of partial derivative of the scalar function (vector) B 
with respect to the elements of vector x. 



All items are published in Russian. 
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